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Abstract 

A sample of 'single-source' events, compatible with the multifragmentation of very 
heavy fused systems, are isolated among well-measured ^^^Gd + ^°-^\] 36 A.MeV 
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reactions by examining the evolution of the kinematics of fragments with Z > 5 
as a function of the dissipated energy and loss of memory of the entrance channel. 
Single-source events are found to be the result of very central collisions. Such central 
collisions may also lead to multiple fragment emission due to the decay of excited 
projectile- and target-like nuclei and so-called 'neck' emission, and for this reason 
the isolation of single-source events is very difficult. Event-selection criteria based on 
centrality of collisions, or on the isotropy of the emitted fragments in each event, are 
found to be inefficient to separate the two mechanisms, unless they take into account 
the redistribution of fragments' kinetic energies into directions perpendicular to the 
beam axis. The selected events are good candidates to look for bulk effects in the 
multifragmentation process. 

Key words: NUCLEAR REACTIONS ""*U(155q^j x), E=36 A.MeV, central 
collisions, selection of "fused" systems. 



1 Introduction 



The study of the behaviour of nuclear matter under extreme conditions of 
temperature and pressure, and the possible associated phase transitions, con- 
stitutes one of the major axes of research in nuclear physics today. The priv- 
ileged experimental tools to conduct this study are heavy ion collisions from 
the sub-relativistic to the ultra-relativistic energy range, coupled with power- 
ful large-acceptance multidetector arrays. With beam energies of a few tens 
to a few hundreds of MeV per nucleon one hopes to explore the environs of 
the 'Hquid-gas' coexistence region, which is predicted to be located at less- 
than-normal densities and at temperatures below the critical temperature for 
the phase transition {Tc ~ 10-16MeV) [1]. In such collisions excited pieces of 
nuclear matter may be formed with excitation energies comparable to or even 
greater than nuclear binding energies. 

In this bombarding energy range, which has been explored at various experi- 
mental facilities around the world over the last decade (e.g. GSI Darmstadt, 
NSCL Michigan, GANIL Caen), an evolution is observed from evaporative pro- 
cesses, in which light particle emission is the principal mode of decay of hot 
nuclei with temperatures < 5MeV, to the so-called multifragmentation regime 
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where large numbers of nuclei with a wide variety of masses are observed in 
the exit channel of individual collisions [2,3]. The goal of multifragmentation 
studies is to link this experimental observation to the properties of the nuclear 
matter phase diagram, which entails disentangling dynamical effects linked to 
the collision phase from decay products of thermodynamically equilibrated 
excited nuclear systems. 

For heavy- ion collisions between 20 and 100 A.MeV incident energy, the binary 
character of almost all reactions has been established experimentally [4-9] ei- 
ther by the identification of surviving projectile-like and target-like fragments 
(PLF and TLF, respectively) or through the observation of reaction prod- 
ucts' characteristics compatible with their origin in the statistical decay of 
two principal primary excited nuclei, the quasi- projectile (QP) and quasi- 
target (QT). It was first thought that such binary reaction mechanisms were 
essentially a continuation of the deep-inelastic collisions well-known at lower 
energies [10,11], accompanied by increasingly important "preequilibrium" or 
"prompt" emissions due to the opening of the phase-space for nucleon-nucleon 
collisions. However, increasingly exclusive experiments with in detector arrays 
have shown that an important fraction of the detected particles and heavier 
fragments originate from the rapidity region between the projectile and tar- 
get [12], and that this contribution is difficult to distinguish from that which 
may be ascribed to the decay of fully equilibrated QP and QT [13,14]. They 
may, depending on the species considered, the beam energy, system size and 
impact parameter, find their origin in e.g. prompt emissions from the overlap 
(participant) zone between the two colliding nuclei [15,16], decay of highly- 
deformed QP/QT [17,18], or the formation and rupture of a necklike structure 
between projectile and target [19-22]. 

Faced with such a highly complex situation for most of the collisions, the 
so-called "single-source" events observed for a small (~ 1%) part of the cross- 
section for collisions of very heavy ions around the Fermi energy [23-25] are of 
great interest. In these events, all of the emitted fragments and particles, apart 
from a small preequilibrium component of light particles, seem to originate in 
the multifragmentation of a single nuclear system containing almost all of the 
available mass and energy of the entrance channel, which simplifies enormously 
the analysis. Such events provide a unique opportunity to study the decay of 
well-defined and very heavy pieces of excited nuclear matter, for which one 
expects that bulk effects, if present (for example the mechanical instability 
associated with spinodal decomposition [26-28]), should play a dominant role. 
Semi-classical transport calculations [29,30] predict these reactions to occur 
for central collisions of heavy nuclei {b < O.Sbmax)- The initial system formed 
should therefore suffer large amounts of compression and heating, and may 
subsequently expand into the low-density coexistence region. 

The method by which a certain set of data is selected depends on the exper- 
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imental conditions, detector systems and the reaction studied. In the present 
work we present in detail the procedure by which single-source events have 
been brought to light and isolated for ^^^Gd + ""*U 36 A.MeV reactions stud- 
ied with INDRA. The method used (based on the angle Oqow signifying the 
polar deviation from the beam direction of the events' principal axis ) has 
already been published elsewhere [23,25], but no comparison has been made 
before between this method and other procedures which are far more com- 
monly used to sort experimental data in the specific context of single-source 
event selection. This is the goal of this paper. In the accompanying paper [31] 
we will use these data in order to probe the role of bulk effects in the multi- 
fragmentation of such heavy systems. 

The experimental set-up, including the detector array specifications and op- 
erating conditions during the experiment, is presented in Sec. 2. The first step 
towards an unbiased global reconstruction of the recorded reactions is to ensure 
that a very large proportion of the emitted products have been measured event 
by event : this is achieved by an initial selection of 'complete events' (Sec. 3.1). 
In order to introduce the comparison between different selection methods, the 
most commonly used global variables are presented in Sec. 3. 2 and 3.3. Then 
the evolution of the reactions from binary collisions to single-source events is 
described in Sec. 4 by classifying them according to ^aow Finally in Sec. 5 we 
show that the most commonly-used selection methods are generally less effi- 
cient for isolating single-source events by comparison with the 6'flow selection. 



2 Experimental details 

The ^'^■''^Gd -|- ""^U system was studied with the Att multidetector INDRA, 
operating at the GANIL accelerator. INDRA, which was described in detail 
in [32,33], can be viewed as an ensemble of 336 telescopes covering ~ 90% of 
the An solid angle. The detection cells are distributed amongst 17 rings cen- 
tred on the beam axis. Low energy identification thresholds and large energy 
ranges were obtained through the design of three-layer telescopes, composed 
of an axial- field ionization chamber operated at 30 mbars of C^Fg, a 300 /xm 
silicon detector and a CsI(Tl) scintillator, thick enough to stop all emitted par- 
ticles, coupled to a phototube. Such a telescope can detect and identify from 
protons between 1 and 200 MeV to uranium ions of 4 GeV. Past 45 °, where 
fast projectile- like fragments are no longer expected, the telescopes comprise 
only two stages, the ionization chamber operated at 20 mbars and the scintil- 
lator. Finally the very forward angles (2 — 3 °) are occupied by Nel02-Nell5 
phoswiches. Charge resolution of one unit was obtained up to Z=64 for the 
fragments identified through the AE — E method in the Si-Csl couple, and up 
to Z=20 when the AE signal is furnished by the ionization chamber (above this 
limit imposed by the energy resolution of the ionization chambers, extrapola- 
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tion of the identification to lieavier fragments is assured by calculations based 
on energy-loss tables, with a resolution of a few charge units). 




Fig. 1. Energy thresholds of the 3-member INDRA modules, full line: charge iden- 
tification thresholds; dashed line: minimum energy to reach the Csl; dotted line: 
detection threshold ( ions deposing 5 Me V in the ionization chamber, but not iden- 
tified). 

Identification thresholds increase from ~ 0.7 A.MeV for Z—3 to ~ 1.7 A.MeV 
for Z—50 (Fig.l). Slow heavy ions at around or just below the threshold energy, 
stopping in the detector directly behind the ionization chamber, are partially 
identified i.e. are assigned a minimum atomic number. In our analysis such 
particles are only included to permit the global characterisation of the least 
violent collisions, when fission fragments of the target nucleus (see Fig. 8(a) 
and accompanying discussion in the text) are present. Energy calibration of 
all detectors ensured an accuracy of within 5% . Great attention was paid 
to pulse-height defect calibration of the Si-dOOfim detectors for fragments 
Z > 15 [34]. For Rings 2-9 residual energy deposited in Csl scintillators 
were derived from energy losses in the (preceding) silicon detectors. Energy 
cahbrations of Csl scintillators were only performed for Rings 10-17 (45 < 
diab < 176°) and were recently improved with a better description of quenching 
which was previously over-estimated because ^-electrons were neglected [35]: 
measured kinetic energies of fragments with Z pa 15-35 detected beyond 45" 
for the ^^^Gd -|- ""*U 36 A.MeV system were noticeably modified as compared 
to previously published data [28] (see accompanying paper [31]). 

A -"^^^Gd beam accelerated by the GANIL cyclotrons was used to bombard 
a 100 ng/crn? U target, enclosed between two 20 ng/crn? carbon layers. In 
order to ensure a negligible rate of multiple interactions in the target, the 
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beam intensity was maintained around 3x10^ pps. The low target thickness 
was chosen to allow slow fragments to escape the target. Due to the high 
detection efficiency of the Csl for 7/e~, a minimum-bias trigger based on the 
multiplicity of fired modules was chosen. Two settings were used: first M>4 
and then M>8 in order to minimise the proportion of triggering events which 
involve the carbon backings. Only runs using the M>8 trigger are analysed 
here. In the off-line analysis, events having a multiplicity of correctly identified 
charged particles inferior to the experimental trigger condition were rejected 
for reasons of coherency (with the exception of Fig. 2). 



3 Event selection methods and procedures 

3. 1 Removal of poorly measured events from the data sample 

As a first step in our event selection procedure, we seek to exclude from the 
data sample poorly-measured events, i.e. ones for which a complete detec- 
tion of all reaction products has not been achieved. This is necessary for two 
reasons. First of all, we want to characterise those rare events for which a 
composite system containing almost all of the mass of the colliding nuclei 
undergoes multifragmentation. Secondly, it is an essential prerequisite for the 
analyses using global variables that we will use later. In the case of the ^^^Gd -|- 
nat^ reaction, a way must also be found to ehminate reactions of the projectile 
with the carbon target backing. 

The degree of confidence with which we can judge the analysed events to be 
'complete' increases with the total detected charge and momentum (Fig. 2). 
Indeed, without making any hypothesis about the physics of the studied re- 
actions, with a 100% -perfect detection one should measure the total system 
charge {Zgyg = Zproj + Zfarg = 64 + 92 = 156) and the beam momentum (Pproj) 
in the exit channel for every event. Fig. 2 shows that experimentally this is far 
from being the case. It should be noted that the 'momentum' used here is cal- 
culated from the product of atomic number Z and velocity component in the 
beam direction Vz, and normalised to the incident (projectile) 'momentum' : 

Ptot ^ T,Zvz 

Pproj ZpfQjVpj'oj 

This procedure is used to compensate for the fact that INDRA does not per- 
mit isotopic identification for fragments Z > 4, and neutrons are not detected. 
Simulations of central collisions for the highly neutron-rich ^"''^Gd + ""*U sys- 
tem including the effect of secondary particle emission show that this 'charge 
momentum' is conserved to around 96% , with a r.m.s. deviation of ~2% . 



6 



Gd + U 36 A.MeV INDRA 




0.5 1 10 20 30 40 

Total momentum P^JP^^^ Total multiplicity 



Fig. 2. Non-equidistant logarithmic contour plots showing event-by- event correlations 
between the total detected charge Z^^x , ^momentum ' Ptot cind multiplicity of charged 
products Nq- The numbered boxes correspond independently to the three classes of 
events discussed in the text. 

The numbered boxes in Fig. 2 correspond to different classes of detected events. 
The design of the INDRA detector was optimised for the study of multifrag- 

mentation, typically high-multiplicity (< 50) events in which reaction prod- 
ucts are emitted over a large angular range. On the other hand, less violent 
collisions are not so well detected due to : 

(1) identification thresholds (see Sec. 2) which exclude heavy, slow-moving 
target residues because of incomplete identification; 

and 

(2) loss of rapid projectile-like fragments at laboratory angles smaller than 
2 °. 

Box 1 events mainly correspond to rather peripheral collisions for which nei- 
ther projectile- (QP) nor target-like (QT) heavy fragments were detected, only 
hght charged particles (LCP, Z = 1,2) and occasionally fission fragments from 
the target, the latter having very low energies which only permit an approx- 
imate identification via a minimum atomic number Zmin (see Sec. 2). Most of 
the events in box 2 have total charge Zfot ~ 70 and total momentum close 
to that of the projectile nuclei. They correspond either to the detection of a 
QP and some hght particles or to close-to-complete detection of ^^^Gd-|--'^^C 
collisions in reverse kinematics. For total detected charge Ztot > 70 reactions 
with the carbon target support are completely excluded. For measured Ztot 
values upwards of 50% of the total charge of the system ^^^Gd -|- ""*U, increas- 
ing 'charge-completeness' basically leads to a global increase of the number 
of detected products without much change to the type of reactions observed 
(except for low-multiplicity events where the principal detected nuclei are a 
quasi- projectile and the target fission fragments). 

Because of a lower overall efficiency for fragment detection compared to LCP, 
and trivially because of fragments' larger atomic numbers, the requirement of 
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highly charge-complete events slightly favours events where a larger proportion 
of the available charge is detected in the form of fragments at the expense of 
light charged particles (see Fig. 3(a) and Figs. 4(e), (f)). This is especially true 



Gd+ " U E/A=36MeV M>8 




Total detected charge Z^^j Total detected charge Z^^j 

Fig. 3. Effects of the total detected charge Ztot of analysed events, (a) mean multi- 
plicities of light charged particles (Nlcp) o.'^'d of fragments with Z > 5 (Nf) as a 
function of Ztot i (b) weighted mean velocity of the Z > 5 fragments (see text) versus 
Ztot- Logarithmic intensity scale representing measured cross section. Points are the 
mean velocities averaged over each Ztot bin. 

for fragments emitted backwards in the cm. frame, because of the increase in 
effective thresholds at large laboratory angles due to the recoil of the cm. . 
Thus the weighted mean velocity of fragments ^miVi/ Y^mi is much larger 
than the ^^^Gd + ""^'U c.m. velocity Vc.m. for small total detected charge, and 
indeed close to that of the projectile, fproj (see Fig.3(b)): we saw above that 
in these events all the backwards-emitted (i.e. target) fragments were missed. 
As Ztot increases the detected fragments' mean velocity decreases as more 
and more backwards-emitted fragments are included. In other words, those 
reaction products which have a low probability of detection are not present 
in 'incomplete' events, whereas in 'complete' events nearly all emitted frag- 
ments must have been detected, including those emitted into experimentally 
unfavourable regions of phase-space, and so the selection of complete events 
appears to favour the detection of those products. 

In order to retain a reasonably-sized data sample of very well measured events, 
we imposed the condition Ztot > 120(~ Q.nZsys) in the off-line analysis. We 
judged this charge-completeness to be sufficient for two reasons. Firstly, for 
Ztot > 120 the weighted mean velocity of detected fragments (see Fig. 3(b)) 
has an average value close to (within 10% ) the c.m. velocity of the collisions 
Vc.m.i which is a necessary condition for all of these fragments to have been 
produced in the disassembly of a fused system containing most of the mass of 
the projectile and target. Angular distributions of fragments and LCP, event 
shapes and orientations (see following section) are identical for this data sam- 
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pie whether calculated in the reaction (^^^Gd + "'^*U) cm. frame or in the 
reconstructed cm. frame. Secondly, if only events with Ztot > O.OZ^j^g are con- 
sidered the only significant difference is that the mean fragment multiplicity 
increases hy ^ 1, while the size of the data sample is drastically reduced. 
Angular distributions and kinematical observables are unchanged. Therefore 
we consider that events with Ztot > 120 are on the average sufficiently well- 
measured to be classified as 'complete events' (box 3 in Fig. 2). The additional 
constraint on the total detected momentum, 0.8 < Ptot/Pproj < 1-1; has little 
effect for these events. 

The measured cross-section (using target thickness and incident ion flux) for 
complete events is 93 mb, to be compared with a calculated reaction cross- 
section value of (Tr = 6.5 barns from systematics [36]. An experimental mea- 
surement of (7r was not possible because of the carbon target support. If this 
selection seems somewhat draconian, let us recall that it is an essential con- 
dition in order to be able to correctly reconstruct the kinematics of events 
where a very large part of the total system has undergone multifragmenta- 
tion. We will show in the following that the accuracy of this reconstruction is 
of paramount importance for the isolation of such events. 



3.2 Impact parameter selectors (IPS) 

How may one isolate a sample of events corresponding to the formation and 
multifragmentation of a single excited nuclear system ? A first answer may be 
to reason in terms of impact parameter. The 'fusion' events that we are looking 
for must correspond to central collisions : in a low energy picture for reasons 
of angular momentum; in a high energy picture to maximise the participant 
zone. 

Events may be classed into impact parameter bins using global variables which 
are supposed to increase (or decrease) monotonically with h and the geomet- 
rical prescription [37] : 



Here $ represents the chosen global variable, in this case reaching its max- 
imum value ^rnax whcu 6—^0. Nf.^ is the total number of recorded events 
(corresponding to the total geometrical cross-section 7r6^^^), and 6e,st is the 
estimated impact parameter for events characterised by the value $ = $i. 
The assumption underlying Eq.( 2) is that one can assign a single impact 
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parameter to each value of the global variable and vice-versa, i.e. that fluctu- 
ations are negligible. 






Fig. 4. (a)-(d) N on- equidistant contour plots showing event-by- event correlations 
between A^lcPj -^C; Etu, Et and Njmf (see text for definitions), (e) Distributions of 
N{ and A'imf for all recorded events (open histogram and triangles, respectively) and 
for complete events (shaded histogram and diamonds, respectively) . (f) Distributions 
of Eti2 for all recorded events (open histogram) and for complete events (shaded 
histogram). An approximate scale in integrated cross-section, beginning with the 
'most central collisions' (largest values of each IPS), is shown in percent for the two 
variables Mmf o-nd -Eti2- 

The most often- used IPS variables are multiplicities (total multiplicity Nq., 
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LCP multiplicity A^lcp) and total transverse energies (of all particles E^, or 
uniquely LCP -E'ti2)- In addition we have defined the multiplicity of 'interme- 
diate mass fragments' (IMF) to be the number of nuclei with 3 < Z < 30 
(A^imf), and the multiplicity of fragments for products with Z > 5. Panels 
(a)-(d) of Fig.4 present event-by-event correlations between these variables. 
One may first note the large spread of values around the mean correlation 
in each case, therefore fiuctuations cannot be neglected, as has often been 
observed for dissipative collisions in and below the Fermi energy range. Each 
IPS will thus share events differently between estimated impact parameter 
bins. On the other hand the effect of fiuctuations on Eq.( 2) is to underes- 
timate the impact parameter for the 10% most central collisions [37-39]. On 
Fig. 4 approximative integrated cross-sections are shown for A^imf and Etu 
which are calculated with respect to the total number of recorded events, irre- 
spective of Zto^. As this number includes reactions between the beam and the 
carbon target backings (see Sec. 3.1) the percentages given are a lower hmit. 
Complete events therefore belong to the > 10% most central ''^^^Gd -|- ""*U 
collisions measured (Fig.4(e) & (f)) and further 'impact parameter selection' 
e.g. multiplicity cuts may well be ineffectual. 

Let us mention in passing the variable Zi,ound used by the Aladin collabora- 
tion [40,12]. It corresponds to the sum of atomic numbers of all fragments with 
Z > 2. In collisions between lOOA.MeV and lAGeV using the Aladin spec- 
trometer Zbound represents the size of the projectile spectator minus evaporated 
hydrogen isotopes. In this energy range the near-geometrical dependance of 
spectator size on impact parameter makes Zi,ound a good IPS. The variable 
Zbound does not behave as an IPS with the complete detection of particles and 
fragments from both target and projectile. In another context, it can be useful 
for model comparisons if it is defined to be the total charge of all fragments 
(e.g. with Z >b) when comparing data with calculations whose principal aim 
is to reproduce fragment partitions but not all the LCP emitted at different 
stages of the reaction (see accompanying paper [31]). For complete events of 
^^^Gd + "°*U the average value of the total charge contained in fragments 
Z >b corresponds to approximately 50% of the total system charge. 

3.3 Global shape variables (GSV) 

A more direct method of discriminating between different reaction mechanisms 
is based on considering how fragments are distributed in the centre of mass mo- 
mentum or velocity space on an event-by-event basis (the 'event shape') [41]. 
We carry out this analysis using fragments with Z > 5 rather than the more 
usual Z > 3 definition because of the very heavy nature of the system (note 
that both definitions are quite arbitrary). LCP are left out of the analysis as 
they result from several different mechanisms not directly related to fragment 
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production (e.g. pre-equilibrium emission, evaporation from hot fragments) 
and also their smaller masses can blur event shapes in velocity space. Then 
the formation of 'fused systems' containing most of the incident nucleons is 
characterised by isotropic emission of fragments in the centre of mass frame, 
assuming for simplicity complete relaxation of the multifragmenting system's 
form i.e. a spherical source, and negligible angular momentum (spin). On the 
other hand the existence of two principal moving sources of emission (e.g. the 
two partners of a deeply-inelastic collision), or of spectator-like fragments sep- 
arating from some 'participant zone' around mid-rapidity, implies a rod-like 
elongated event shape. The variables used to exploit this difference in event 
shape come from high-energy particle physics, where pattern recognition in 
momentum distributions has long been a major concern [42]. 

From the cartesian components of fragment {Z > 5) momenta in the centre 
of mass one may construct the tensor [43], 

Z>5 



whose eigenvectors and eigenvalues may be interpreted in terms of an ellipsoid 
in momentum space, and which give information on both the event shape 
(sphericity, coplanarity etc.) and its orientation with respect to the beam (see 
Fig. 5, upper panels). If {Ai, A2, A3} are the eigenvalues of Eq.(3), normalised 
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Fig. 5. Values of global shape variables ( GSV) for different event-shapes (represented 
by an ellipse) and orientations with respect to the beam axis, 'b' and 'P. A. ' refer to 
variables calculated with respect to the beam axis or the principal axis of the event, 
respectively. It should be noted that spherical event shapes (central panel) can only 
be achieved for very high (infinite) multiplicities. 

to their sum and labelled according to increasing size, A3 > A2 > Ai, then 
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sphericity -5=^(1 -As) (4) 
a/3 

coplanarity C = -— (A2 — Ai) (5) 

The 'flow axis' corresponding to the largest eigenvalue A3, or major axis of 
the eUipsoid, defines the principal axis of the event, whose deviation from 
the beam is referred to by the polar angle ^aow This angle varies between 
0° and a maximum value of 90° because, by convention, it is taken as the 
angle between the forward lobe of the ellipsoid and the beam. For events 
which have lost all memory of the entrance channel kinematics (e.g. single- 
source events with negligible angular momentum) this angle is isotropically 
distributed, whatever the fragment multiplicity. In the case of a single source 
with large angular momentum all angles are populated but directions close to 
the beam are favoured [39], due to the preferential emission of fragments in 
the reaction plane. 

Independently of the tensor one may construct other shape-dependent vari- 
ables, for example : the second moment of Fox and Wolfram, H2 [44] 



^ E^IpiI IP2I (3cos2^,e^-l), (7) 



H(Q) ^ 2 
1 
2 



m-Y.\\pA \P2\ (8) 



1,2 



where the sums are over pairs of fragments, pi,P2 are their cm. momenta, and 
Orel is their relative angle; the momentum isotropy ratio, Riso [45], 

^iso = ^(Eb"'J/Ekl|) (9) 



where the sums are over all fragments, and p±,p\\ represent perpendicular and 
parallel projections of p; and an energy isotropy ratio such as Eiso [46] , 

E,,,^1-^(^EJY.E^ (10) 
£^± = E^iSin2^i (11) 

i 

where the sums are again over fragments, Ei is the cm. kinetic energy of 
the i^^ fragment, and 6i its polar angle. An alternative choice, with a slightly 
different definition, is the variable ERAT used by the FOPI collaboration [47]. 
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By construction H2 is independent of the choice of basis axes in the cm. mo- 
mentum space (it only depends on relative angles between pairs of fragments), 
whatever the event shape. The values of E^so and Ri^o for spherical events are 
also independent of the choice of axes. However, spherical events require infi- 
nite multiplicities of fragments even if the latter are isotropically emitted, and 
so are never observed [48]. This is because the event-shape deduced from a 
small number of isotropically-emitted fragments is only a poor approximation 
to the 'true' spherical distribution, due to the extremely restricted sampling 
of the available momentum space event by event. As Fig. 5 shows, because the 
event shape is always non-spherical E^so and i?iso have to be calculated with 
respect to the principal axis of the event in order to ensure non-ambiguous 
shape determination. When calculated with respect to the beam axis their val- 
ues depend not only on the event shape but also its orientation in the velocity 
space, ^flow In the following we will implicitly suppose that isotropy ratios are 
calculated unambiguously, unless explicitly stated otherwise. 

The effectiveness of event-shape discrimination is weakened by the fact that 
we cannot unambiguously define what values of the GSV should correspond 
to the formation and decay of fused systems because for small multiplicities 
the event-shape depends strongly on the number of emitted fragments. On 
the other hand, an isotropic distribution of flow angles for a set of events 
is an unambiguous signal which is independent of the fragment multiplicity 
that fragments have lost all memory of the entrance channel. This fact will be 
exploited in the following section in order to select single-source events. 

Fig. 6 presents some experimental correlations of these global shape variables 
for complete events (an event-shape analysis has no sense unless performed for 
complete events). Fig. 6(a) shows the correlation between a global shape vari- 
able, /2iso, and an impact parameter selector, -E'ti2- Two regions are apparent: 
one corresponding to very elongated events and quite large impact parame- 
ters {R\so < 0.25, £'ti2 < 300 MeV); the other covering a very broad range of 
both event shapes and transverse energies. The former set of events are the 
least dissipative collisions which are sufficiently violent in order to be well- 
measured i.e. which are not excluded by the requirement of complete events. 
We will see in the next section that the principal products of these reactions 
are residues of the quasi-projectile and of the target fission fragments which 
are well-separated in velocity space. As for the remaining events, they corre- 
spond to strongly dissipative collisions and it should be noted that the event 
shape and impact parameter observables show no strong correlation. We will 
show the existence of 'fusion' i.e. single-source events for a fraction of these 
collisions in the following section. 

Figs. 6(b) and (c) show the experimental correlations which exist between the 
three GSV. All of them well separate the two classes of events remarked in 
Fig. 6 (a), especially H^. Esso shows a greater variation than H2 for the most 
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Fig. 6. Event-by- event correlations, for complete events, between: (a) event-shape 
(Riso) and impact parameter selector Eti2; (b) E\so and H2; (c) R\so and 
Eiso(isotropy ratios calculated with respect to the principal axis of each event); (d) 
R\so calculated with respect to the beam axis and with respect to the principal axis 
of each event. 

dissipative collisions, which may mean that it has a greater sensitivity to vari- 
ations among very compact event-shapes. On the other hand the two isotropy 
ratios, calculated from fragment momenta or kinetic energies, are very strongly 
correlated and may be considered to give equivalent information on event 
shapes. 

Finally in Fig. 6(d) we present the correlation between the values of the isotropy 
ratio i?iso calculated with respect to the beam axis ('ambiguous' shape vari- 
able) or the principal axis of the event ('unambiguous' shape variable). For 
the most elongated event shapes (small -Riso) the two methods of calculating 
Riso give very similar values, because the least dissipative events have mean 
6'flow angles of ^10° (see Fig. 7), therefore the principal axis and the beam 
axis are almost aligned. The majority of the more compact events (principal 
axis isotropy ratio > 0.25) also have axis-independent isotropy ratios (crest 
following a line "beam axis -Riso=principal axis i?iso"), and this is once again 
because 6'flow remains quite small for these events (Zone 2 of Fig. 7). 

However there are also events for whom the isotropy ratio (i.e. the apparent 
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event shape) depends strongly on the axis chosen, for example appearing as 
"spherical" (/?iso = 1) with respect to the beam axis but non-spherical {Riso < 
1) with respect to the principal axis of the event. Note also the existence of 
events for whom the isotropy ratio with respect to the beam axis it!iso > 
which implies an event shape elongated perpendicular to the beam direction. 
In both of these cases the 'ambiguous' apparent event shape is due to the 
increase of ^^aow, which modifies the projection of the fragment momenta on 
to the beam axis. For a given unambiguous Riso approaching 1, the truly 
increasingly spherically-symmetric events populate 'beam axis event shapes' 
which vary less and less with ^aow, as one would expect. However we do not 
observe an 'island' of truly spherical events with -R^l""* = ^fso^' = 1- Once 
again this is because for the fragment multiplicities in play, a spherical event 
shape can not be attained due to finite number effects, even if all fragments 
are emitted isotropically. 



4 Event classification using the 'Wilczynski diagram': isolation of 
single-source events 

In order to classify ^^^Gd + ""*U 36 A.MeV coUisions we use a method which 
was first employed for the analysis of Pb-|-Au 29 A.MeV reactions studied with 
the multidetector Nautilus [23]. The correlation between total measured cm. 
kinetic energy of detected charged products (TKE) and the principal direction 
of fragment 'flow' (^how) when plotted for each event (see Fig. 7) resembles the 
Wilczynski diagram [49] well-known at bombarding energies below 20 A.MeV 
. TKE provides a dissipation scale for the collisions, decreasing as more and 
more excitation energy is deposited in the system. For complete events we 
may write 

TKE = Ecm. + Q ~ ^ Eneutron ~ E "^7 i^'^) 

where Ecm., Q, ^Eneutron ^^d T^E^ are, respectively, the available centre- 
of-mass energy, the mass balance of the reaction and the total neutron and 
gamma ray kinetic energies. The increase of E* with increasingly dissipative 
collisions implies higher average multiplicities and energies of neutrons, light 
charged particles, and fragments : the latter account for an increasingly neg- 
ative mass balance. 

Just as at low energy (for heavy systems), a dog-legged correlation appears 
with a crest running from forward-peaked, slightly dissipative collisions (large 
TKE) to highly-damped reactions with little or no memory of the entrance 
channel (large 6'flow)- We have defined four classes of events depending on their 
position in the 'Wilczynski diagram', Fig.7. In order to examine the topology 
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Fig. 7. 'Wilczynski diagram' for complete events: logarithmic intensity scale repre- 
senting measured cross-section as a function of total measured cm. kinetic energy 
(as a fraction of the available centre of mass energy ) and 'flow ' angle or direction of 
the event principal axis, 0flow The four zones indicated are used to classify complete 
events (see text). For events in Zones 2 to 4 the mean value of TKE is indicated 
(points) for each 6'flow bin. 

of the fragment emission in velocity space in each class of events, we present 
the charge-velocity correlations for each Zone, Fig.8. Here, as a function of 
velocity along the beam axis or the principal axis of each event, the average 
charge density of fragments is reported, which we define as the average charge 
of all the fragments with Z > 5 which fall in a given bin in parallel velocity, 
normalised to the width of the bins. This is a simplified version of the tool 
first presented in [50]. The distributions of Fig.8 therefore give an insight into 
the repartition in velocity space of the charge bound in fragments. 

In Zone 1 are found the least dissipative collisions which are sufficiently violent 
in order to be well measured. The direction of the principal axis for these events 
remains close to the grazing angle (6'^^™' ~ 9-^° [36]), and the principal reaction 
products are residues of the quasi-projectile {v ~ 4.5 cm/ns) and of the target 
fission fragments {v ~ —2 cm/ns and -4.5 cm/ns, respectively). Let us note 
in passing that although Fig.8 (a) shows clearly that the majority of charge 
bound in fragments is concentrated in projectile- and target-like products, the 
velocity space between projectile and target (the so-called mid-rapidity region) 
is not empty but populated by a few intermediate mass fragments which may 
result from the decay of a 'neck' of matter between the two principal partners 
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Fig. 8. Fragment charge-velocity correlations for the four classes of events corre- 
sponding to the zones defined in the Wilczynski diagram (Fig. 7). The mean charge 
density (see text) is plotted as a function of the velocity parallel to the beam axis 
(thick lines) or the principal axis, P. A., of each event derived from the tensor Eq.3 
(shaded histograms). 



of the reaction, as has been observed for mid-central colhsions of other heavy 
systems around the Fermi bombarding energy (see for example [19-21,13]). 

The most dissipative collisions (Zones 2 to 4) populate all 6'flow angles, although 
in majority events are still concentrated at small angles, ^aow < 30° (Zone 2 in 
Fig. 7). The cos ^flow-distribution becomes fiat (isotropic) for large flow angles 
(Fig.9(a)). The repartition of the fragments' charge is more homogeneous than 
in Zone 1: specifically, the mid-rapidity region is fully populated. However, 
in Zone 2 (Fig.8(b)) there are two asymmetric 'bumps' in the distribution at 
forward and backward velocities, strongly reminiscent of the entrance channel. 
Although not strictly correct, we will refer to this type of fragment emission 
topology as 'binary' events in the following. In Zones 3 and 4 (Fig.8(c) and (d)), 
due to the increasing deviation of the principal axis from the beam direction, 
6'flow, the charge density is not distributed in the same way according to one or 
the other axis. Thus, with respect to the beam direction, the two components' 
relative velocity becomes smaller in Zone 3 than in Zone 2 while in Zone 4 only 
one component peaked at the cm. velocity is seen; on the other hand, in the 
velocity-space along the principal axis, two components are present whatever 
the Zone and their relative velocity decreases only very slightly while the 
distribution becomes symmetric with regard to the cm. velocity. 
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This result is difficult to interpret. The distributions with respect to the beam 
axis suggest that as Oqow increases fragments are produced in the break-up of 
increasingly relaxed projectile- and target-like primary fragments, giving way 
to fusion-like reactions for Oq^w > 70°. However, this would be in contradiction 
with the fact that the mean TKE of events in Zones 3 and 4 is constant, and 
only slightly lower than in Zone 2 (Fig. 7, points). The energy dissipated in the 
collisions is practically constant for Zones 2 to 4, which is, on the other hand, 
consistent with the velocity distributions with respect to the principal axis of 
each event. One can then understand the apparently spurious evolution of the 
two components' relative velocity along the beam direction to be due to the 
effect on the projections on to the beam axis of increasing ^flow- 
Does this then mean that fragments in Zone 4 result from the break-up of 
projectile- and target-like primary fragments, moving apart almost perpen- 
dicularly to the beam ? This is not so clear, as in fact both distributions of 
Fig. 8(d) are compatible with the break up of a fused system in to a small num- 
ber of fragments. This is because event-shapes for low-multiplicity reactions 
are never spherically symmetric (see Sec. 3. 3), therefore one can always define 
event-by-event an oriented principal axis. In the case of the isotropic break- 
up of a fused system the principal axis has no physical significance and can 
take all directions, therefore angular distributions with respect to the beam 
are isotropic. On the other hand the principal axis lies by definition in the 
direction of maximum elongation of each event shape, therefore with respect 
to this axis fragments seem to be preferentially emitted in the forward and 
backward directions, with equal probability [39]. Thus Fig.8(d) does not per- 
mit to conclude on the origin of the fragments, and we must examine their 
characteristics in greater detail. 




Fig. 9. For the most dissipative collisions (Zones 2^4 of Fig. 8): (a) 'flow' angle 
distribution; (b) mean sphericity of the events as a function of the fragment multi- 
plicity Nf; (c) distribution of relative angles Orel between pairs of fragments in Zone 
4 events; (d)-(g) evolution of fragment- fragment relative angle distributions with 
O&ow, ecLch distribution has been divided by the distribution for Zone 4 events. 
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Mean event-shape sphericity as a function of fragment multiphcity is shown 
in Fig. 9(b) for the three Zones in ^aqw corresponding to the most dissipative 
colhsions (Zones 2-4). One should first of all note that for each class of events 
the increase of sphericity with fragment multiplicity is expected from the finite 
number effects we discussed above. We see that in general event sphericity 
increases with 'fiow' angle, and that for each fragment multiplicity the most 
compact events are found in Zone 4 (^^aow > 70°). This is a strong signal of the 
evolution of fragment kinematics towards single-source emission, as : (i) the 
event sphericity and 6'flow are a priori independent observables; and (ii) only 
comparisons of event shapes for constant multiplicity avoid distortions due to 
the aforementioned finite number effects. 

One may examine in more detail this evolution by looking at the distribu- 
tions of relative angles {9rei) between pairs of fragments emitted in the same 
event. Fragments emitted by the same source have an isotropic distribution 
of 9rei angles in the rest frame of the emitter, except for small 9rei which are 
suppressed by Coulomb repulsion between nascent fragments in the case of 
rapid successive emissions. Fig. 9(c) shows exactly this type of distribution for 
fragment-fragment relative angles in the Zone 4 events, which are therefore 
compatible with the fast break-up of a single source. Figures 9(d)-(g) show the 
evolution of the relative angle distributions with the flow angle of the events. 
In each case, the distribution has been divided by that of the Zone 4 events 
(Fig.9(c)), in order to highlight any differences between the two event samples 
(a value of 1 for all angles means that the considered distribution is identical to 
that for Zone 4 events) . The fact that small and large relative angles between 
fragment pairs are favoured compared to Zone 4 in events with ^aqw < 10° 
(Fig.9(d)) reveals the dominance of emission from two distinct sources mov- 
ing apart in the cm. frame in these events: large relative angles are populated 
by fragment pairs coming from different sources, while the repopulation of 
small relative angles is due to kinematical focusing in the direction of the 
source velocity of fragments born of the same source. The decreasing relative 
population of large and small angles when increasing the flow angle of the 
events reflects the diminishing importance of collisions leading to projectile- 
and target-like primary fragments. It should be noted that relative angle dis- 
tributions such as that shown in Fig. 9(f) can also be obtained if fragments are 
emitted from a spherical source having very large angular momentum [39] . 

To summarise the results presented in this section, we have shown that the 
'Wilczyhski diagram' (Fig. 7) allows to sort well-measured reactions in to 
classes of events ('Zones') according to the total measured kinetic energy in 
the cm. frame TKE and the flow angle ^aow We have presented the evolution 
of these reactions from the least violent collisions, where flssion of the target 
nucleus is observed, to highly dissipative 'binary' collisions, in which many 
fragments are produced but they retain a strong kinematical memory of the 
projectile-target asymmetry. This evolution concerns the flrst three Zones of 



20 



the 'Wilczynski diagram'. In Zone 4 are found events for which : (i) no memory 
of the entrance channel remains; (ii) no projectile- or target-like fragments are 
observed; (iii) mean event shapes are the most compact of all the very dissi- 
pative collisions; (iv) fragment kinematics are consistent with rapid emission 
from a single source; (v) moreover this single source must contain a very large 
proportion of the available mass and energy i.e. preequilibrium emission is 
limited (otherwise for an asymmetric system in direct kinematics one would 
expect the 'fused' system to recoil with a velocity less than that of the system 
centre of mass, cf. Fig. 8(d)) [51]. 

The 'Wilczynski diagram' therefore permits the isolation of a sample of single- 
source events for ^aqw > 70". Let us recall that for single-source events the flow 
angle takes all values, therefore single-source events must also be present for 
^flow < 70°. However as we have seen, events with small 6'flow are dominated 
by a 'binary' fragment emission topology, and keep a strong memory of the 
entrance channel. Only at large flow angles do single-source events become 
dominant and separable using a 6'flow cut. Of course, this selection method 
can only function if the flow angle retains a memory of the direction of the 
primary projectile- and target-like fragments coming from highly dissipative 
binary reactions, and this direction must remain close to the beam axis. It is 
possible that this may only occur for heavy systems, for which binary highly 
dissipative collisions are strongly focused around the grazing angle due to 
Coulomb effects [23,25]. If binary collisions of light systems lead to large flow 
angles even with low cross-sections then single-source events, when present, 
may not be revealed by a ^flow-cut. Other more refined selection methods must 
be apphed [52] in this case. 

To conclude this section we will give an estimate for the cross-section associ- 
ated with the formation and multifragmentation of fused systems for ^^^Gd -|- 
""*U 36 A.MeV reactions. The measured cross-section for events in Zone 4 is 
2.6 mb, which is ~3% of the total cross-section for complete events. Now we 
must estimate the fraction of all single-source events represented by this sam- 
ple taking into account (i) the 6'flow selection and (ii) the selection of complete 
events. 

If all the single-source events have an isotropic fiow angle distribution the 
fraction of these events that we have selected is given by cos 70° = 0.34. This 
is the case assuming the formation of spherical fused systems with negligible 
angular momentum. In a less ideal scenario we should allow for deformed 
and/or turning fused systems formed at finite impact parameters. The cos ^flow- 
distribution in this case would be forward-peaked [53] and the fraction of 
single-source events retained by the ^flow cut would therefore be less than 34% . 
This fraction must however be much greater than 3% (fraction of all complete 
events found in Zone 4) because Zone 2 and Zone 3 events are clearly not 
dominated by single-source events. Let us suppose as a conservative estimate 
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that Zone 4 events correspond to 20% of all single-source events. Wc will use 
this as an upper estimate for the total cross-section. As a lower estimate we 
will suppose that the 6'flow-distribution is indeed isotropic and that the Zone 4 
event sample is polluted by 'binary' events. Given the characteristics of these 
events presented above we think that an upper estimate for this pollution is 
around 10% . 

As far as the selection of complete events is concerned, these single-source 
multifragmentation events belong to the class of very violent collisions for 
which the efficiency of the INDRA detector array is optimal. However the 
geometrical efficiency of 90% is perhaps not attained due to the very high 
multiplicities of these events (double hits). The distribution of Z^ot for Zone 4 
events is a monotonously decreasing function of Zi^t (it is of course truncated 
at Ztot = 120 because of our complete events selection) and can be fitted 
by the tail of a gaussian function of mean value < Zt^t >— 115 and width 
a — 13. We estimate the fraction of all single-source events retained by the 
complete-event selection to be the fraction of the area under this gaussian 
curve contained between the limits Ztot = 120 and Ztot = 156 which we found 
from numerical integration to be 37% . 

With these hypotheses we can estimate an upper limit for the total single- 
source cross-section (assuming a deformed spinning source, assuming the Zone 
4 sample to be free of 'pollution' and correcting for event completeness) to be 

<Zrce - 2.6mb X -i^ X -i- = 35mb (13) 



while a lower limit (assuming a pollution of Zone 4 by binary events, assuming 
only spherical zero-spin fused systems and correcting for event completeness) 
is given by 



<Zrce = 90% X 2.6mb X ^ X ^ = 19mb (14) 



Let us remark that if one makes the extreme assumption that all of the most 
dissipative complete events are compatible with the multifragmentation of 
deformed single sources [54] (3% instead of 20% in Eq. (13)) then the upper 
limit for the cross section associated with the multifragmentation of fused 
systems would be estimated to be 234 mb (Ri3% of the total reaction cross 
section) . 



22 



5 Can single-source events be isolated using other selection meth- 
ods ? 



We have shown how a sample of single-source events may be isolated based 
on the degree of rotation of the event ellipsoid, 9fiow This sample represents 
at most 34% of all the single-source events present among the most dissipa- 
tive events. One might wonder if another approach, perhaps based on the more 
commonly-used tools of impact parameter selection or event-shape discrimina- 
tion discussed in Sec. 3. 2 and Sec. 3. 3, would not be better adapted to isolate a 
larger data sample e.g. by including the small-6'flow single-source events ? With 
this goal in mind we will now look in more detail at the behaviour of IPS and 
GSV with regard to the 'Wilczyhski' event classification scheme. 
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Fig. 10. (a) 'Wilczynski diagram' (Fig. 7) for complete events, showing the definition 
of the event position variable, fl. Points correspond toQ = 0, .1, .2, . . .. (b) Evolution 
of total (Nc), LCP (Ni^cp), Z > 5 (N^) and 3 < Z < 30 (Nimf) multiplicities as 
a function of event position in the 'Wilczynski diagram', (c) Evolution of the total 
transverse energy of all charged products, Et, and of LCP, -E'ti2- (d) Evolution of 
the event shape according to H2, Eiso o,nd Riso(the two isotropy ratios are calculated 
according to the principal axis of each event). Points correspond to mean values, 
vertical bars to standard deviations. 

Fig. 10 presents the evolution of the mean values of the IPS and GSV variables 
presented in Sees. 3.2 and 3.3 as a function of event position in the 'Wilczynski 
diagram'. In order to do this we defined an observable, Q, which varies as a 
function of TKE and 6'flow in such a way that it increases monotonously along 
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the crest of the experimental correlation observed in Fig. 7 and Fig. 10(a). fl is 
therefore a global variable related to the position of events in the 'Wilczyhski 
plot'. Its value is given by 



for 8" < ^flow < 14° and 0.48 < TKE/^Je.m. < 0.8, covering values < fl < 0.5, 
and by 



for 14° < ^flow < 30° and 0.48 < TKE/Ec.m. < 0.54 (0.5 < Q < 0.6) and 
for 30° < ^flow < 90° and 0.45 < TKE/Ecm. < 0.53 (0.6 < fl < 1). Outside 
of these limits (represented by rectangular boxes in Fig. 10 (a)) Vl is undefined 
and the corresponding event is not included in the calculated mean. Roughly 
speaking, events from Zone 1 of Fig. 7 have < Q < 0.25; from Zone 2, 
0.25 < O < 0.6; from Zone 3, 0.6 < O < 0.85; and from Zone 4, 0.85 < ^ < 1. 

It is clear that the majority of the evolution of the IPS and GSV variables 
takes place for fl < 0.4 (Fig.lO(b)-(d)) i.e. Zones 1 and 2 of the 'Wilczynski 

diagram', and this evolution is therefore strongly correlated to the dissipation 
of the available kinetic energy Ec.m.- We have shown that these events retain 
a strong memory of the entrance channel (Fig. 8(a), (b) and Fig.9(d),(e)). IMF 
and LCP multiplicities (Fig. 10(b)) as well as the total transverse energies 
(Fig. 10(c)) all increase as the collisions become more dissipative, while event 
shapes become more and more compact (Fig. 10(d)) reflecting the decreasing 
relative velocity between primary target- and projectile-like fragments in the 
exit channel and the increasingly isotropic emission of fragments in the cm. 
frame. 

If wc now look to the most dissipative collisions {Q > 0.4: Zones 2-4 of Fig. 7) 
the mean values of IPS variables show a slight evolution, much smaller than the 
standard deviations of their distributions (vertical bars in Fig.lO(b)-(d)). Only 
Et continues to increase significantly for Q > 0.5 because increasing ^gow means 
increasing fragment transverse energies; however Et is constant for D, > 0.8. 
The mean Z > 5 fragment multiplicity A^f is the same whatever the position 
of the event in the horizontal branch of the 'Wilczyhski diagram'. It should 
be noted that IMF and LCP multiphcities show an identical saturation with 
dissipated energy, unhke what was observed for the ^^^Xe-|-^°^Bi 28 A.MeV 
reaction [55]. 

IPS are therefore quite insensitive to the evolution of fragment kinematics 
from 'binary' collisions to single-source events that was shown in Sec. 4, and 
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they seem to be most strongly correlated with energy dissipation. Let us point 
out that a simple calculation for heavy systems of the excitation energy per 
nucleon for a fully-damped pure binary collision (no mass exchange) 



E* = E, 



'cm. 



1 ApAt 
2Ap + At 



(17) 



(where Vrei is given by the Viola systematics [56]) or for fusion 



E* = E,,^, + A{Ap, Zp) + A(A, Zt) - A{Ap + At, + Zt) 



(18) 



(where A{A, Z) is the mass excess of nucleus using the extrapolation 
of [57] in the case of the fused system) gives very similar results. This may 
explain the lack of significantly different IPS behaviour for single-source events. 

Event-shapes on the other hand continue to evolve towards more compact 
forms with increasing Ofiow (Fig. 10(d), Q > 0.6) as in Fig.9(b), but here we can 
sec that, as for IPS, distributions about the mean values are comparatively 
wide. It should be noted that although here we mix events with different 
fragment multiplicities, the mean Z > 5 multiplicity is the same for all Q > 0.4 
(Fig. 10(b)) and so the evolution of the GSV truly reflects a change in event 
shapes, i.e. the evolution of the fragment kinematics with increasing ^aow This 
evolution is complete for Oqow ^ 70° {ft > 0.85): it is for this reason that we 
choose to define our sample of single-source events (Zone 4 of Fig. 7) as having 
6'flow > 70° even though the distribution of relative angles between fragments 
for events with 60° < 6'flow < 70° is identical to that for ^aow > 70° (see 



5.1 Selecting the most central collisions 

Let us try to define IPS cuts corresponding to very central collisions in order to 
isolate single-source events. In Fig. 11 the IPS- variable distributions for Zone 
4 events (shaded histograms) are compared to those for all recorded events 
(including ^^^Gd -|- ^^C collisions, lines) and for the most dissipative events 
(Zones 2 to 4 of Fig. 7, hatched histograms). All of these variables increase 
monotonously with decreasing impact parameter, i.c the largest values corre- 
spond to the most central collisions (see Sec. 3. 2). Fig. 11 shows that Zone 4 
events correspond to central collisions (large IPS values) but even though they 
have mean values of IPS which are higher than for Zones 2 and 3 (Fig.lO(b,c)), 
the distribution is so wide that the same values are explored by Zone 4 events 
as by all of the most dissipative events. This observation can be interpreted 
in terms of a large overlap of impact parameters leading to single-source or to 
'binary' exit channels, with single-source events merely a subset of the most 



Fig.9(g)). 
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Fig. 11. Comparison of IPS distributions for all recorded events (histogram), the 
most dissipative events (Zones 2-4 of Fig. 7, hatched histogram) and Zone 4 events 
(shaded histogram). All histograms are normalised with respect to the total number 
of recorded events. 

central collisions. It is clear that no cut on an IPS variable can be defined 
which would isolate some or all of Zone 4 events without mixing them with 
the other Zones. The only exception is Et (Fig. 11(e)) for which, because large 
^fiow angles imply large fragment transverse energies, Zone 4 events are the 
only ones to explore the largest Et values. 

One may then wonder if high-multiplicity or high transverse energy cuts can 
be used in order to isolate a sample of single-source events different to that 
constituted by Zone 4 of the 'Wilczyhski diagram'. Such cuts have been applied 
to the most dissipative events i.e. Zones 2 to 4 of Fig. 7 (we have shown that 
IPS are capable of distinguishing between these events and those of Zone 1, 
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Fig. 6(a)) and the results are presented in Fig.l2. The cuts are : A^c ^ 37, 
A^LCP > 29, > 1070MeV, Etu > 570MeV. These values are slightly higher 
than the mean values observed for Zone 4 events (cf. Fig. 11), and they were 
chosen in part so that the number of events retained by each cut was the same 
within 10%. 

Flow angle distributions for the events selected with A^"c, iVLCP and Etu cuts 
(Fig. 12(a), (d),(j), shaded histograms) are forward-peaked and appear to be of 
the same form as the total distribution (thick lines) except that the proportion 
of events with the smallest ^aow is shghtly reduced. We see that only ~10% 
of the Zone 4 event sample is retained by these IPS cuts. It should be noted 
that for ^^^Au+-'^^''Au 35 A.MeV reactions (very similar in both total mass and 
available energy to the studied reactions), a much smaller total multiplicity 
cut (A^c > 24) was found to be sufficient in order to isolate a sample of highly 
spherical events with an isotropic flow angle distribution corresponding to 
10% of the total measured reaction cross-section [24]. Our results show that 
such a cut is completely insufficient for reactions of an otherwise identical 
asymmetric reaction. The cut applied to total transverse energies (Fig. 12(g)) 
on the other hand greatly reduces the proportion of events with cos ^aqw > 0.5 
and the resulting ^aq^ distribution is nearly isotropic. The proportion of Zone 
4 events retained by this cut is larger (~30% ) than for the other IPS as is to 
be expected from comparison of Fig. 11(e) with Fig.ll(a)-(d) and (f). 

If any of these cuts isolates a sample of single-source events, the characteristics 
of the events should be independent of the flow angle. Therefore we separate 
the events selected with each IPS cut into two lots, those with small flow 
angles (^^aow < 70°) and those with large flow angles (^aow > 70°). It should 
be recalled that the second lot is a subset of the Zone 4 event sample. 

Distributions of relative angles between fragments are shown in Figs. 12(b), 
(e), (h) and (k). In each case, the distribution has been divided by that of 
the Zone 4 events (Fig.9(c)), in order to highlight any differences between 
the two event samples (a value of 1 for all angles means that the considered 
distribution is identical to that for Zone 4 events). For the multiplicity and 
LCP total transverse energy cuts we see that in the sub-sample of events with 
^flow < 70°, large and small fragment-fragment angles are favoured compared 
to Zone 4 events (thick lines) suggesting fragment emission from two sources. 
This means that the Onow distributions for events selected by these cuts are 
forward-peaked because they mix single-source and 'binary' events. On the 
other hand the events selected by the £'t-cut show fragment-fragment relative 
angle distributions which are identical to Zone 4 events independently of the 
subsequent division according to flow angle. 

The average charge density of fragments in the different sets of events are 
presented as a function of their velocity along the principal axis in Figs. 12(c), 
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Fig. 12. Results of applying IPS cuts defined in the text to the most dissipative colli- 
sions (Zones 2-4 of Fig. 7). (a),(d),(g),(j): flow angle distribution of the events se- 
lected by the IPS cut (shaded histograms) compared to the total distribution for Zones 
2-4 (thick lines). (b),(e),(h),(k): distribution of fragment-fragment relative angles 
for IPS-selected events divided bin-by-bin by the distribution obtained for Zone 4 
events (Fig. 9(c)) i.e. R{cos9rei) = 1 for identical distributions. (c),(f),(i), (I): frag- 
ment charge-velocity correlations for IPS-selected events with respect to the principal 
axis of each event. In the latter two cases, the events selected by IPS cut are also 
separated into those with 0fiow < 70° (thick lines) and those with 0flow ^ 70° (shaded 
histograms) . 

(f), (i) and (1). The 6'flow < 70° subset of the events selected by cuts on A^c, 
A^LCP and i?ti2 shows a memory of the entrance channel through the presence of 
larger fragments at velocities v\\ < 0. Note also that these cuts seem to favour, 
among the events with 6'flow > 70°, events having a larger mean charge density 
for v\\ > (heavier fragments emitted in the forward direction). The E^- 
selected events present very homogeneous, symmetric and 6'flow-independent 
charge-velocity correlations, although a slight increase of the mean charge 
density is perhaps visible at backwards velocities for 6'flow < 70° (Fig.l2(i), 
thick line). 

To sum up this section on the selection of the most central collisions, we 
have shown that the sample of single-source events found in Zone 4 of the 
'Wilczyhski diagram' is a subset of the most dissipative events as far as IPS are 
concerned, and are therefore probably a subset of the most central collisions. 
Attempts to isolate single-source events with an IPS cut supposed to select 
the most central collisions only result in mixed samples containing single- 
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source and 'binary' events. The exception to this rule is the total transverse 
energy which selects a sample of single-source events with an almost isotropic 
flow angle distribution, albeit with some pollution by 'binary' events being in 
evidence for the smallest Oqow 



5.2 Selecting the most compact event shapes 

Let us now consider the selection of the most spherical (or compact) event 
shapes, by defining cuts using global shape variables. Fig. 13 compares the 
event shapes of Zone 4 events with those of all the most dissipativc events. 
For all GSV (except Eiso calculated with respect to the beam axis) most of the 
values populated by all of the most dissipative events can also be associated 
with Zone 4 events. In all cases it can be seen that the Zone 4 single-source 
event sample explores a very wide range of event shapes. One may then worry 
that cuts made to restrict to only the most compact event shapes will select a 
very particular subset of single source events (see below). As in the IPS case, it 
is clear that no cut on a GSV variable can be defined which would isolate all or 
some of the Zone 4 events without mixing them with the other Zones, except 
for the two 'ambiguous' isotropy ratios. For these GSV Zone 4 events are the 
only ones to explore values signifying a large proportion of fragment emission 
in directions perpendicular to the beam (see Sec3.3, Fig.6(d)), because these 
events have the largest 6'fiow angles. 

We have defined GSV cuts supposed to select the most spherical event shapes, 
and in each case the size of the event sample was the same (within 10% ) as 
those previously selected with IPS cuts. The cuts were: H2 < 0.06, -Eiso < 0.21 
and Riso > 0.74 for the 'unambiguous' shape variables, |£'iso| < 0.5 and 0.93 < 
i?iso < 1-07 for the isotropy ratios calculated with respect to the beam axis. 
First, let us consider the results for the 3 'unambiguous' GSV (Fig. 14). 

Selecting the most compact event shapes docs not change greatly the form 
of the flow angle distributions (Fig. 14(a), (d), (g)). Events with small Oqow 
are still favoured in the selected event samples, although the Eiso cut does 
slightly flatten the cos ^aow distribution. The fragment-fragment relative angle 
distributions, on the other hand, arc independent of ^^aow for each GSV cut 
(Fig. 14(b), (e), (h)). Moreover the suppression of small and large 6rei relative to 
the distribution for Zone 4 events is the opposite of what one would expect for 
samples mixing single-source and 'binary' events. In fact, these distributions 
are directly related to the GSV cuts as can be readily seen by inspection of 
the H2 definition, Eq.8. As H2 depends on cos'^Oreh imposing small values of 
H2 favours events where the majority of fragment pairs have |cos^^re/| ~ 0. 
Although the link between 9rei and the isotropy ratios is not so trivial to 
demonstrate analytically, it is very probable that the same kind of effect is in 
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Fig. 13. Comparison of GSV distributions for the most dissipative events (hatched 
histogram) and single-source events (19flow ^ 70°, shaded histogram). E[so(beam) and 
Riso(beam) are the values calculated with respect to the beam axis (cf. EisJ^ and Riso'' 
of Fig. 5), Eiso and Riso are calculated with respect to the event principal axis (cf. 
Eiso^'^' and Riso^"^' of Fig. 5). All histograms are normalised with respect to the 
total number of complete events. 

operation also. 

This is not the only non-trivial effect that cuts on 'unambiguous' shape vari- 
ables can have on the event topology. The charge-velocity correlations show 
that the Riso and (to a lesser degree) H2 cuts favour events in which the heav- 
iest fragments have small velocities parallel to the principal axis. This may 
be a way of 'compacting' the event in momentum space. The charge- velocity 
correlations for i^iso-selected events on the other hand are very homogeneous, 
as one would naively expect for highly compact event shapes. 
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cos 6jj^^ cos 6^|,| V// principal axis (cm/ns) 

Fig. 14. Results of applying GSV cuts defined in the text to the most dissipa- 
tive collisions (Zones 2-4 of Fig. 7). (a),(d),(g),(j): flow angle distribution of the 
events selected by the GSV cut (shaded histograms) compared to the total distribu- 
tion for Zones 2-4 (thick lines). (b),(e),(h),(k): distribution of fragment-fragment 
relative angles for GSV-selected events divided bin-by-bin by the distribution ob- 
tained for Zone 4 events (Fig. 9(c)) i.e. R{cos9rei) = 1 for identical distributions. 
(c),(f),(i), (I): fragment charge-velocity correlations for GSV-selected events with re- 
spect to the principal axis of each event. In the latter two cases, the events selected 
by GSV cut are also separated into those with O^ow < 70° (thick lines) and those 
with 0fiow > 70° (shaded histograms). 

Finally let us see what is the effect of applying to the most dissipative collisions 
cuts designed to select the most compact events using the two 'ambiguous' 
isotropy ratios, E^^o and i?iso calculated with respect to the beam axis. The 
resulting event samples do indeed appear to be good samples of single-source 
events: they have fragment-fragment relative angle distributions which are 
identical to that of Zone 4 events whatever value takes 6'flow (Fig-15(b),(e)) and 
charge-velocity correlations show an homogeneous, symmetric distribution of 
the fragments in velocity space (Fig. 15(c), (f)). These cuts therefore avoid the 
peculiar effects on event topology that are seen with the 'unambiguous' GSV. 
However Fig. 15 (a), (d) shows that the flow angle distribution of the selected 
events is unphysical, having been distorted by the correlation between the 
'ambiguous' isotropy ratios and 6'flow Small ^flow events are strongly suppressed 
or even excluded (in the case of -Eiso) by the 'ambiguous' GSV cuts. 

To summarise this section on selecting the most compact events, we have 
shown that the single-source event sample isolated in Zone 4 of the 'Wilczyhski 
diagram' explores a very wide range of event shapes, therefore GSV cuts 
can only select subsets of single-source events with very particular fragment 
emission topologies. Specifically, the 'unambiguous' GSV select events with 
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Fig. 15. Results of applying 'ambiguous' GSV cuts defined in the text to the most dis- 
sipative collisions (Zones 2^4 of Fig. 7). (a), (d): flow angle distribution of the events 
selected by the GSV cut (shaded histograms) compared to the total distribution for 
Zones 2-4 (thick lines). (b),(e): distribution of fragment- fragment relative angles for 
GSV-selected events divided bin-by-bin by the distribution obtained for single-source 
events (Fig. 9(c)) i.e. R{cos9rei) = 1 for identical distributions. (c),(f): fragment 
charge-velocity correlations for GSV-selected events with respect to the principal 
axis of each event. In the latter two cases, the events selected by GSV cut are also 
separated into those with 0flow < 70" (thick lines) and those with 6'flow > 70° (shaded 
histograms) 

forward-peaked flow angle distributions, in which pairs of fragments are pref- 
erentially detected at 90° and the heaviest fragments are at rest in the centre 
of mass frame. On the other hand the isotropy ratios when calculated with 
respect to the beam axis permit the extraction of small-^Aow single-source 
events which are similar to the Zone 4 event sample. However most of the 
single-source events with small 6'flow (cos^aow < 0.3) remain irrecoverable due 
to the strong correlation between 6'flow and these variables. 



6 Conclusions 



We have shown for ^^^Gd + ""^U 36 A.MeV reactions measured with the IN- 
DRA 47r detector array how a sample of single-source events may be isolated 
which correspond to the formation and multifragmentation of very heavy fused 
systems, comprising the majority of the entrance channel nucleons. These 
events are selected from among the most well-characterised reactions (at least 
80% of the total charge and momentum were measured) using a condition 
on the flow angle ^aqw between the principal axis of the event ellipsoid, con- 
structed from fragment {Z > 5) kinetic energies, and the beam direction. One 
expects events with large 6'flow to show little memory of the colliding nuclei 
because of the forward-focused differential cross-section for deeply inelastic 
collisions of heavy nuclei. We have shown that a class of events compatible 
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with the fast break-up of a compact single source is dominant at large flow 
angles, ^aow > 70". 

This sample of single-source events is associated with the largest mean mul- 
tiplicities and total transverse energies of charged reaction products, and on 

average the most compact event shapes. However the single-source events are 
only a subset of the most central/most compact events, therefore isolating a 
sample of them using impact parameter selectors (IPS) or global shape vari- 
ables (GSV) is not as trivial as it may first appear. We found that the only 
variables capable of isolating an event sample with the same characteristics 
as events with Oq^w ^ 70°, namely the total transverse energy of charged 
reaction products, Riso and i^iso the isotropy ratios of momentum and kinetic 
energy flow, respectively, calculated with respect to the beam axis, are all cor- 
related with the transverse kinetic energy of the detected fragments, and this 
is also true for 6'fiow Therefore using the present methods it is not possible to 
obtain an unbiased ^gow distribution (which could give important information 
on the relaxation in form of the multifragmenting sources, angular momentum, 
etc. [53,54]) for the totality of the single-source events in the data sample. 

The ^flow > 70" cut permits the study of a sample of very heavy multifrag- 
menting systems, for which bulk properties may play a decisive role. This work 
is presented in the accompanying paper [31]. 
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